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We present the -D4 x Z 2 flavor symmetry, which is different from the previous work 
by Grimus and Lavoura. Our model reduces to the standard model in the low energy 
and there is no FCNC at the tree level. Putting the experimental data, parameters are 
fixed, and then the implication of our model is discussed. The condition to realize the 
tri-bimaximal mixing is presented. The possibility for stringy realization of our model is 
OO ! also discussed. 
O 

1. Introduction 

X 

It is the important task to find an origin of the observed hierarchies in masses and flavor 
mixing for quarks and leptons. Neutrino experimental data provide us an important clue 
for this task. Especially, recent experiments of the neutrino oscillation go into the new 
phase of precise determination of mixing angles and mass squared differences [TH2] . Those 
indicate the tri-bimaximal mixing for three flavors in the lepton sector [3]. Therefore, it 
is necessary to find a natural model that leads to this mixing pattern with good accuracy. 

Flavor symmetries, in particular non-Abelian discrete flavor symmetries, are interesting 
ideas to realize realistic patterns of mass matrices. Actually, several types of models with 
non-Abelian discrete flavor symmetries have been proposed [3J. Furthermore, non-Abelian 
discrete flavor symmetries can be realized in the simple geometrical understanding of 
superstring theory [ofo] as well as extra dimensional models. The _D 4 symmetry can appear 
typically in heterotic string models on factorizable orbifolds including the Z 2 orbifold. 
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Indeed, several semi-realistic models with D4 flavor symmetries have been constructed 
in Ref. [5117] and in those models three families correspond to a singlet and a doublet 
under the flavor symmetry. Therefore, taking _D 4 symmetry as the flavor symmetry of 
quarks and leptons, these mass spectra and the flavor mixing matrix should be carefully 
examined to establish the realistic model of quarks and leptons [8 9J. 

The -D4 flavor symmetry was at first proposed for the neutrino mass matrix by Grimus 
and Lavoura [8j . In this model, the atmospheric neutrino mixing is maximal while the solar 
neutrino mixing is arbitrary They introduced three electroweak Higgs doublets together 
with two neutral singlets in the scalar sector to reproduce the large flavor mixing angles. 
Then, the tree level flavor changing neutral scalar vertices do not vanish. Moreover, when 
we consider supersymmetric extension of this -D4 flavor model, such a supersymmetric 
model would have three pairs of up and down Higgs fields. That would violate the gauge 
coupling unification, which is one of important aspects of the minimal supersymmetric 
standard model, unless one introduces extra colored supermultiplets0 

In this paper, we propose alternative -D4 flavor model with one Higgs doublet, which 
reduces to the standard model in the low energy. There is no tree level flavor changing 
neutral current (FCNC) in our model. The higher dimensional operators provide the 
charged lepton and neutrino masses. Putting the experimental data, our parameters are 
fixed, and then the implication of our model is discussed. 

The paper is organized as follows: we present the framework of the -D4 model in Sec. 
2, and discuss the neutrino masses, flavor mixing angles and Higgs potential in Sec. 3. 
In Sec. 4, the numerical results are discussed. Section 5 is devoted to the summary and 
discussion. 

2. D4 flavor symmetry and Yukawa couplings 

We present the framework of our D 4 flavor model. The D4 symmetry has five irreducible 

representations, that is, a doublet 2 and four singlets, 1 ++ , 1_ + and 1 , where 1 ++ 

is a trivial singlet and the others are non-trivial singlets. Their products are decomposed 

as 

2x2 = l ++ + l + _ + l_+ + l__, l ab xl cd = l ef , (1) 

where a,b,c,d = ±, e = ac and / = bd. Here, the left-handed lepton doublets are denoted 
as I a (a = e, /1, r) and the right-handed charged leptons and right-handed neutrinos are 
denoted as e#, /j,r,tr, N e , N^, N T , respectively. The first family leptons are assigned to D4 
trivial singlets, while second and third family ones are to D± doublets. The electroweak 
Higgs doublet h is a -D4 trivial singlet. We summarize charges of flavor symmetry in Table 
[TJ where new gauge singlet scalar fields x, Xii X2 are introduced and additional Z 2 
charges are assigned for leptons and scalars. 

2.1. Charged lepton mass matrix 

We write down the Yukawa interactions, which are invariant under the gauge group of 
the standard model and the flavor symmetry D 4 x Z 2 , by using the multiplication rule of 



7 We would study a supersymmetric D4 model in a separate paper [ID] . 



D4, Savor symmetry 



3 







h 


X X-+ {xi,X2) 


z 2 


1++ 2 1 ++ 2 1 _|__|_ 2 

+ + - + + 


1++ 

+ 


1++ 1-+ 2 
+ 



Table 1 

-D4 and charges given for leptons and scalars. 



D A in Eq. © 
-A = 



?/ e e/ e /lX + 2/t(/^ + + Vr(.-plp + Tl T )hx- + 

+Vere{ l nXl + lrX2)hx + y'^ilpXl ~ lrX2)hX- + 
+yre(p>Xl + rX2)khx + VreifiXl + T^O^X-H- 

+^r(Ml + TX 2 )(V%1 + l rX2)hx + vlriflXl ~ f X2)(^Xl ~ k^O^X 

+y c ^r(fiX2 + fxi)(lpX2 + Wxi)hx + ytr(fiX2 - rxi)(l»X2 - kxi)hx 

+y' a ur(fiXi + r% 2 )(LXi - KX2)hx-+ + y'iAfiXi - r% 2 )(Lxi + l T X2)hx~+ 



flT' 

+y'%{UX2 + rxi)(^X2 - lrXi)hX-+ + y'ir(UX2 - fXi)(^X2 + Wxi)hx-+ 
+ --- + h.c, (2) 

where x an d X-+ denote x/A and X-+/A respectively, and A is the cutoff scale. The scale 
A is taken to be the Planck one in our numerical study. The ellipsis in Eq. (jSJ) denotes 
higher order contributions but they are negligibly small in our considerations. 
We take the vacuum expectation values of scalar fields as follows: 



(h)=v, (( X i,X2)) = (V 1 ,V 2 ), ( X )=V a , ( X -+) = V b 



(3) 



where v = 174GeV and others are taken to be D4 symmetry breaking scale. After 
spontaneous symmetry breaking, the mass matrix of charged lepton becomes 



Mi 



y e a a 



(ye^a - yL a b)a (y e ^ a + JJ ' OL h ) 



a 



{y^oi a - y'^abja y T a a - y' T a b {y^o-a + y'^b)® 
(y^ e aa a + y'a b )a (y MT a a - y' a b )a 2 y T a a + y' T a b 



(4) 



where a a = V a /A and a b = V b /A and we assume the vacuum alignment in the D4 doublet 
scalar field, V\ = V2, so that, ((xi> X2)) = (V, V). The parameter a is defined as a = V/A. 
This vacuum alignment is important for the masses and mixings in the neutrino sector. 
Since the value of a is sufficiently small as discussed later, the charged lepton mass matrix 
can be approximately regarded as diagonal. The masses of charged leptons are given by 



m e = y e a a v, 



m, 



y T a a v - y' T a b v, m T = y T a a v + y' T a b v. 



(5) 



We need the fine-tuning to obtain the difference between the masses of the muon and the 
tau, m^/m r <C 1, as discussed in Ref. [8]. 

2.2. Neutrino mass matrix 

Let us consider the neutrino sector. We can write down the possible Dirac mass terms 
up to the dimension five operators by the same prescription as the charged lepton sector, 
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-£ D = yiNJ e h + y 2 (N^l^ + N T l T )h 

+ynN e (l fl Xi + lrX2)~h + VzxiNpXi + N T Xi)leh 
+yUN„X\ + N T X2)(lf,Xi + lrX2)h 

+VU^Xl - N T X2)(hdl - lrX2)h 
+023 (ftfo + fXl)(i M X2 + lrXl)h 
+1/23(^X2 - fXl)(l^X2 - lrXl)h 

+ ■■■ + h.c. , 

where h = ir 2 h* . The Majorana mass terms are given as 
C R = M 1 N^C~ 1 N e + M 2 (N^C~ 1 N I1 + N^C~ X N T ) 
+y a NjC- 1 (N llX i + N T X2) 
+y a b (Nlxi + N? X 2)C-\N»Xi + K X 2)/A 
+V b b (N^xi - Nj X 2)C~\N^xi - n tX2 )/a 
+V c b (N^X2 + NjxJC-'iN^ + N tX i)/A 

+y d b {Nlx2 - Njxi)C-\N^ X 2 - N tX i)/A 
+ ■■■ + h.c. 

Then the neutrino mass matrices of Dirac and Majorana Mr are given by 



(6) 



M, 



v 



yi yva 

y%\OL y 2 
y 21 a y 32 a 2 



yi2d 

y23® 2 
y2 



M, 



Mi y a Aa y a Aa 
1 y a Aa M 2 y b Aa 2 
y a Aa y b Aa 2 M 2 



(7) 



(8) 



Similarly to the case of charged leptons, the ellipses in Eqs. (j§ 
higher order contributions but they are negligible. 

The neutrino mass matrix is given by the see-saw mechanism, 

M v = M D M^ X M T D . 

The neutrino mass matrix has the following structure, 



M,, 



and d2J) correspond to 



(9) 




(10) 



where 



.4 
C 



y\M\ 



M X M 2 - 2a 2 A 2 M 2 yl ' 

y 2 {M x M 2 -y 2 a a 2 A 2 ) 

M l M 2 -2a 2 y 2 a A 2 M 2 



B 
D 



yiy 2 y a aAM 2 



Aft M 2 -2a 2 y 2 a A 2 M 2 

y&WA 2 

M 1 M 2 -2a 2 y 2 a A 2 M 2 



In these expressions, higher order terms are neglected under the assumption of 
MiMf > a 4 A 3 , M 2 > a 2 M 1 , M 2 > a 2 A, A > M 1 . 



(11) 



(12) 



D4 flavor symmetry 



5 



These assumptions are justified by the numerical analyses as discussed later. The neutrino 
mass matrix is diagonalized by the following mixing matrix, 

/ c s \ 

V = -s/y/2 c/V2 l/y/2 , (13) 
V s/V2 c/V2 -1/V2 J 

where c = cos 9\ 2 and s = sin 6* 12 and 9\ 2 corresponds to the solar mixing angle [8] . Then 
the neutrino mass matrix Eq. (1101) is represented by the solar mixing and neutrino mass 
eigenvalues (i — 1 ~ 3) such as M v = V A diag(m 1 , m 2 , m 3 )V T , which is 




c 2 mi + s 2 m2 — cs(rrii — m 2 )/\/2 — cs(mi — m 2 )/v / 2 

— cs(mi — m 2 )/y / 2 (s 2 mi + c 2 m 2 + m 3 )/2 (s 2 mi + c 2 m 2 — m 3 )/2 
— cs(mi — m 2 )/y / 2 (s 2 m! + c 2 m 2 — m 3 )/2 (s 2 mi + c 2 m 2 + m 3 )/2 



(14) 



and we have the relations, 



Av 2 = c 2 mi + s 2 m 2 , 5f 2 = — ^L(mi — m 2 ), 

V2 

1 1 
6f = — (s mi + c m 2 + m.3), .Dtj = — [s mi + c m 2 — 777,3). (15) 

For neutrino masses, we find 

777i+ 777 2 = (A + C + D) V 2 , 

V^2 2 

777i — 777 2 = Bv , 

CS 

m 3 = Cv 2 - Dtj 2 . (16) 
Then, the mass squared differences and the solar mixing angle are expressed by 

Am 2 tm = -1(a + C + D-^b\ 7/ + (C - D) V, 

Ami, = (A + C + D) ^Bv\ 

cs 

™t29 12 = -L-(C-A + D). (17) 

2.3. Potential analysis 

Here, we analyze the scalar potential and discuss the assumption of vacuum alignment, 
V\ = V 2 . The relevant scalar potential of (x, X-+> Xii X2) is given by 

-c v = -nix 2 - vix 2 - + - l4(Xi + xl) (18) 

+A1X 4 + \ 2 xt + + A 3a (x? + xl) 2 + hbixl ~ xl) 2 + hcXlxl 

+Ai2X 2 X 2 + + Ai 3 X 2 (Xi + xl) + A23X 2 +(Xi + xl) + Ai23XX-+(Xi ~ xl)- 
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The minimum conditions are 

° C = 2V a + 2\ ± V 2 + A 12 V^ + A 13 (n 2 + V 2 2 ) + A 123 ^(Vf - K 2 ) ) = 0. 



a 



2V b + 2A 2 V, 2 + A 12 \/ 2 + A 23 (Vf + if) + A 123 ^(lf - V 2 ) ) = 0, 



j£ = 2V 1 [-^ + 2\ 3a (V 2 + V 2 ) + 2A 3b (F 2 - V 2 ) + A 3c lf 
+Ai 3 K 2 + A 23 lf + A 123 KH ) = 0, 
2V 2 ( - /i 2 + 2A 3a (V? + V 2 2 ) - 2X 3b (V 2 - V 2 ) + A 3c lf 
+Ai 3 K 2 + A 23 H 2 - A123KH ) = 0. (19) 
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Since there are sixteen parameters (^1,2,3, Ai )2i3a)36)3c , Ai 2 ,i3,23,i23, K,6,i, 2 ) while there are 
four equations, these equations can be solved. For this analysis, the following relation is 
important, 

(4A 36 - A 3c )(lf - V 2 ) + 2A 123 KH = 0, (20) 

which is derived from dC v /dxi = and dC v /d\2 = 0. To align the vacuum of V\ = V 2 , 
one requires Ai 23 = 0, which is an assumption in our model. We may impose additional 
symmetry to realize A i23 = 0. Inserting A i23 = 0, we have 



* a 



(2A 2 A 3 — 2A 23 )^i 2 + (2Ai 3 A 23 — Ai 2 A 3 )/z 2 + (2Ai 2 A 23 — 4A 2 Ai 3 )^t| 



4AiA 2 A 3 + 4Ai 2 Ai 3 A 23 — 4AiA 23 — 4A 2 A 13 — A 3 A 12 
2 (2A i3 A 23 — Ai 2 A 3 )/i 2 + (2AiA 3 — 2A 2 3 )/i 2 + (2A i2 Ai 3 — 4AiA 23 )/i 2 



4AiA 2 A 3 + 4Ai 2 Ai 3 A 23 — 4AiA 23 — 4A 2 A 13 — A 3 A 12 
(Ai 2 A 23 — 2A 2 Ai 3 )^t 2 + (Ai 2 Ai 3 — 2AiA 23 )^ 2 + (4AiA 2 — A 2 2 )/x| 

l2 ~7T7T2 _ \ _\2 



4AiA 2 A 3 + 4Ai 2 Ai 3 A 23 — 4AiA 23 — 4A 2 A 2 3 — A 3 A 2 2 



where A 3 = 4A 3a + A 3c . It is found that we can take V a ~ V b , which is necessary to obtain 
muon and tau masses by adjusting parameters. 



3. Numerical discussion 



Let us discuss our numerical results. We define the following two dimensionless param- 
eters, which are the ratios of M 2 and a A to Mi, respectively, 
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By using these parameters and Eq. (TTTT) , the mass squared differences and the solar mixing 
angle are rewritten as 



Am 



-{y\r 2 + y\r + V^ymy gkr / 'csf / '4 + y\{r - 2y 2 a k 2 ) 2 v 4 
~ <r 2 -2y 2 a k 2 r) 2 ' M 2 ' 



Am, 



2 -V2y l y 2 y a kr(ylr 2 + y\r) v 4 



r 2 - 2ylk 2 r) 2 cs M 2 
cot 2^ = ^~ y V (23) 

The neutrino masses are given as 

1 y\r 2 + y\r + \f2yiy 2 y a kr / cs t; 2 
mi ~ 2 ' r 2 -2y 2 a k 2 r * Mi' 

1 y2 r 2 + y 2 r . _ ^/2 yi y 2 y a kr / CS V 2 

m2 ~ 2 ' r 2 - 2yjk 2 r * M x ' 

_ y 2 2 (r - 2y 2 A: 2 ) W 2 

3 ~~ — 2 o 2; 2 — XT' ^ ' 

— 2y£k r M 

When we put the best fit values of mass squared differences and the solar mixing angle 
as Aml tm = 2.4 x 10~ 3 eV 2 , Am 2 ol = 7.6 x 10" 5 eV 2 , and sin 2 12 = 0.32 % we have typical 
values of parameters in this model, 

r = 0.12, k = -0.80, M l = 4.9 x 10 15 GeV, (25) 

where we take all Yukawa couplings as yi = y 2 = y a = yi = 1. By taking the cutoff scale 
A as the Planck scale 2.43 x 10 18 GeV, we find 

\a\ = 1.6 x 10~ 3 . (26) 

Therefore, the assumption to regard the diagonal mass matrix (j3J) are justified. The 
assumption of Eq. ffT2|) turns to 



\r\ > \a\ 2 , \r\ 2 > \ak\ \r\ > \ak\, \k\ > |a| , (27) 

which are also justified by the result in Eq. (1251) . The neutrino masses are given as 

mi ~ 3.4meV, m 2 ~ — 9.4meV, m 3 ~ 49meV , (28) 

which indicate the normal mass hierarchy. 

In the above numerical results, we have assumed all Yukawa couplings to be 1. Now let 
us consider how much the above results change by varying Yukawa couplings. Following 
the above results, we assume that |r| <C 1 and ky a = 0(1) for yi,y 2 = Then, we 

approximate Eq. f[2"3l as 

a 2 vi v * a 2 V2yiyl v A y 2 . . 

~ pns* ■ ~ -ijjPCMj • cot2 "' 2 ~-57w' (29) 



3 Note that either k or y a can be small because only their product ky a appears in Eq. ([23]) . 
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Figure 1. Semilogarithmic plots for a versus y\. 



Hence, the parameters k, r, and M\ are obtained as 

k ~ ^l^tan2^ 12 ~0.9xA, r „ /^k (^Y ~ 0.2 x (^Y , 

M x ~ 2v/2V f CQt A 32 ^ 12 ^ V2 ^ 2 ~3x 10 15 x GeV, (30) 
V. cs Amf ol y 

which leads to 

Mik y x y 2 

a = — - — ~ 0.001 x . (31) 

A y a 

The ratio 2/2/2/1 must be of 0(1) in order that the above approximation is valid, i.e. 
y a k = 0(1). Thus, values of k, r and Mi are of the same order as those in Eq. (|25|) . 
However, the value of a would change its order in some region even if we vary yi, yi and 
y a by 0(1), because a depends basically on a cube of 0(1) parameters, i.e. 2 3 ~ 10 
and 0.5 3 ~ 0.1. Let us investigate this behavior numerically. We use Eq. ( 1231) and vary 
yi, y 2 , and y a in the range of 0.5 — 1.4 and taking account for the 3a error-bar of input 
experimental data Am^ tm , Am s 2 ol) and 6*12. We show the random plots of a versus y\ in 
Fig. [U It is found that the value of a is predicted around 10~ 4 — 10~ 2 . The dependences 
of the value of a on other Yukawa couplings such as y 2 and y a are similar to the case of 
y\. Thus we obtain small a as long as Yukawa couplings are of 0(1). 

4. Summary and Discussion 

We have presented the x Z 2 flavor symmetry, which is different from the previous 
work by Grimus and Lavoura. Our model has one Higgs doublet although the neutrino 
mass matrix has the same structure as the one in the model by Grimus and Lavoura. Our 
model reduces to the standard model in the low energy and there is no FCNC at the tree 
level. 
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In order to realize the tri-bimaximal mixing, the condition cot 29u = 7^ must be 
satisfied. Then, we have the condition y\r — y\ = y\y%y„k- Taking Yukawa couplings to 
be order one, this condition turns to simple one r ~ k + 1, which is easily realized by 
adjusting parameters in our model. 

It would be interesting to study supersymmetric extension of our model. In such a su- 
persymmetric D4 model, we would have a specific pattern of superpartner mass matrices. 
We would study it in a separate paper [TO] . 

Finally, we comment on the possibility for stringy realization of our model. The D4 
flavor symmetry can be derived e.g. from heterotic string models on factorizable orbifolds 
including the Z 2 orbifold like Z 2 X Z M orbifolds [ofB] . Indeed, several semi-realistic models 
have been constructed with three families 07], where three families consist of D4 trivial 
singlets and doublets. From this viewpoint, our D4 flavor structure would be natural. 
However, such orbifold models include only D4 trivial singlets and doublets, but not 
-D4 non-trivial singlets as fundamental states. The D4 non-trivial singlet X-+ plays an 
important role in our model. We need to assume that X-+ is a composite scalar of 

doublets, in order to obtain \ h from the Z 2 orbifold. Another possibility would be 

factorizable heterotic orbifold models including the Z^ orbifold like Z^ x Zu orbifolds, 
because such orbifold models can lead to the D4 flavor symmetry, where non-trivial singlets 
as well as trivial singlets and doublets can appear as fundamental modes [5J. Thus, it 
would be interesting to consider the realization of our D 4 model from Z4 orbifold models. 
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